CERN-TH/97-258 
PSI-PR-97-30 
hep-ph/97T0521 



S3 



in 



o 



Non-factorizable photonic corrections to 
c + c — ► WW — ► 4 fermions 



oo 

On 

^ ■ A. Denner 



Paul-Scherrer-Institut, Wiirenlingen und Villigen 
CH-5232 Villigen PSI, Switzerland 



S. DlTTMAIER 

m _ 

>■ \ Theory Division, CERN 

fSj . CH-1211 Geneva 23, Switzerland 

in 



M. Roth 



OS ■ Paul-Scherrer-Institut, Wiirenlingen und Villigen 

CH-5232 Villigen PSI, Switzerland 

JpH. and 

Q-i| Institut filr Theoretische Physik, ETH-Hdnggerberg 

^ CH-8093 Zurich, Switzerland 
> 

x. 

Abstract: 

We study the non-factorizable corrections to W-pair-mediated four-fermion production 
in e + e~ annihilation in double-pole approximation. We show how these corrections can be 
combined with the known corrections to the production and the decay of on-shell W-boson 
pairs, and how the full off-shell Coulomb singularity is included. Moreover, we find that 
the actual form of the real non-factorizable corrections depends on the parametrization 
of phase space, more precisely, on the definition of the invariant masses of the resonant 
W bosons. For the usual parametrization the full analytical results for the non-factorizable 
corrections are presented. Our analytical and numerical results for the non-factorizable 
corrections agree with a recent calculation, which was found to differ from a previous one. 
The detailed numerical discussion covers the invariant-mass distribution, various angular 
distributions, and the lepton-energy distribution for leptonic final states. 
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1 Introduction 



After the successful experimental investigation of the Z boson at LEP1, at present the 
W boson is studied via its pair production at the upgraded collider LEP2. Since W bosons 
decay very rapidly, the actual reaction observed is e + e~ — > 4fermions. The virtual 
W bosons can be reconstructed by selecting those events in which the invariant masses of 
the final-state fermion pairs are close to the W-boson mass. The expected integrated lumi- 
nosity of LEP2 is 500 pb -1 JTJ. Accordingly, one expects to produce about 8000 W pairs, 
and the typical experimental accuracy is of the order of one to a few per cent. For the 
measurement of the W-boson mass the expected accuracy is even below the per-mille level 

This experimental accuracy should be matched or better exceeded by the precision of 
the theoretical predictions. The calculation of the cross section for e + e~ — > 4 fermions with 
an accuracy of 1% or better is a non-trivial task. In the LEP2 energy region, the lowest- 
order cross section is dominated by the diagrams that involve two resonant W bosons. 
All other lowest-order diagrams are typically suppressed by a factor Tw/Mw, but may be 
enhanced in certain regions of phase space. All these contributions are required at the 1% 
level (see Ref. and references therein). At the level of 0(a) corrections, on the other 
hand, it is sufficient to take into account those contributions that are enhanced by two 
resonant W propagators (doubly-resonant contributions). These can be split in two types: 
factorizable and non-factorizable corrections Q. The factorizable corrections are those 
that correspond either to the production of the W-boson pair or to the decay of one of the 
W bosons. The corrections to on-shell W-pair production || and decay || are available. 
The non-factorizable corrections are those in which the production subprocess and the 
decay subprocesses are not independent; typically an additional particle is exchanged 
between two of the subprocesses. Simple power-counting reveals that among the non- 
factorizable corrections only those Feynman diagrams contribute to the doubly-resonant 
contributions in which a soft photon is exchanged between the subprocesses. 

In this paper we explicitly calculate these non-factorizable photonic corrections. It has 
been shown that they vanish in inclusive quantities, i.e. if the invariant masses of both 
W bosons are integrated out |7| . However, for non-inclusive quantities these corrections 
do not vanish in general. The non-factorizable photonic corrections have already been 
investigated by two groups. Melnikov and Yakovlev || have given the analytical results 
only in an implicit form and restrict the numerical evaluation to a special phase-space con- 
figuration. Beenakker, Berends and Chapovsky have provided both the complete formulae 
and an adequate numerical evaluation 0, |I0"|| , but do not find agreement with all results 
of Ref. || . For this reason, it is worth-while to present the results of a third independent 
calculation. 

We start out by discussing the definition of the virtual and real photonic non-factoriz- 
able corrections in double-pole approximation in detail. Since only soft photons are rele- 
vant in double-pole approximation, the virtual non-factorizable correction is just a factor 
to the lowest-order cross section. For the corresponding real correction, the situation is 
similar, but in addition an integration over the photon momentum has to be performed. 
This requires a specification of the phase-space parametrization, which includes, in par- 
ticular, the invariant masses of the W bosons. Usually these are defined via the invariant 
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masses of the respective final-state fermion pairs and are chosen as independent variables 
||, |§ . Experimentally, however, the invariant mass of a W boson is identified with the 
invariant mass of the associated jet pair that necessarily includes soft and collinear pho- 
tons. Therefore, the influence of the choice for the invariant masses of the W bosons on the 
non-factorizable corrections should be investigated in order to provide sound predictions 
for physical situations. 

Besides the non-factorizable doubly-resonant corrections, the most important effect of 
the instability of the W bosons is the modification of the Coulomb singularity. Since the 
off-shell Coulomb singularity results from a scalar integral that also contributes to the 
doubly-resonant non-factorizable corrections, it seems to be natural to approximate this 
integral in such a way that both effects are simultaneously included. This requires going 
beyond the strict double-pole approximation. 

The paper is organized as follows: in Sect. |] we discuss the basis of the double-pole 
approximation, with particular emphasis on the real corrections, and formulate a gen- 
eral gauge-independent definition of the non-factorizable doubly-resonant corrections. In 
Sect. ^ we sketch the calculation of the scalar integrals and describe the evaluation of the 
5-point functions in detail. Section |] summarizes our analytic results for the correction 
factor; in particular, the cancellations between virtual and real photonic corrections are 
inspected, and the off-shell Coulomb singularity is embedded in the approach. Section [5| 
is devoted to the discussion of numerical results. The appendices provide more informa- 
tion on the evaluation of the real bremsstrahlung integrals and the explicit results for the 
required scalar integrals. 

2 Definition of the approximation 

2.1 Conventions and notations 

In this paper we discuss corrections to the process 



e + ( P+ )+e-(p_) -> W + (A; + )+W-(M -> h{k x ) + f 2 {k 2 ) + f 3 (k 3 ) + h(h). (2.1) 



The relative charges of the fermions fi are represented by Qi with i — 1, . . . , 4. The masses 
of the external fermions, mf = k\ and m\ = p±, are neglected, except where this would 
lead to mass singularities. The momenta of the intermediate W bosons are defined by 



their complex mass squared and their respective invariant masses are denoted by 



k. 



+ 



h + k 2 , 



k 






(2.3) 



respectively, and we introduce the variables 



a: 



+ 




M 2 , 



K. 



M 2 . 



(2.4) 
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Furthermore, we define the following kinematical invariants 

s = (p+ + p_) 2 = {k + + k_) 2 , 
S{j [ki ~\~ kj J j 

Sijk = {ki + kj + k k ) 2 , i, j,k = 1,2,3,4, (2.5) 
which obey the relations 

S = k\ + k 2 _ + S13 + S U + S 23 + S 2 4, Sl2 = ^+ , «34 = 

Sjjfc = + Sjfc + Sjfc, i, j, k = 1, 2, 3, 4. (2-6) 
2.2 Doubly-resonant virtual corrections 

The aim of the present paper is to evaluate the non-factorizable corrections to the process 
Q2.1D in double-pole approximation (DPA). The DPA takes into account only the leading 
terms in an expansion around the poles originating from the two resonant W propagators. 

In DPA, the lowest-order matrix element for the process fl2.1|) factorizes into the matrix 
element for the production of the two on-shell W bosons, A4|orn^ W+W (P+>P-> 
the (transverse parts of the) propagators of these bosons, and the matrix elements for the 
decays of these on-shell bosons, Ai^ OT ^^ 1 ^' 2 {k + , k±, k 2 ) and M^ OT ^^^ i {k^, fc 3 , k 4 ): 

KA ST^ - /W Born - /w Born - /w Born / 7 \ 

-MBorn = 2^ TTl? • 

A + ,A_ K + K - 

The sum runs over the physical polarizations \± of the W ± bosons. 

The higher-order corrections to (|2.1|) can be separated into factorizable and non- 
factorizable contributions {§. In the factorizable contributions the production of two 
W bosons and their subsequent decays are independent. The corresponding Feynman 
diagrams can be split into three parts by cutting only the two W-boson lines. The cor- 
responding matrix element factorizes in the same way as the lowest-order matrix element 



The non-factorizable corrections comprise all those contributions in which W-pair pro- 
duction and/or the subsequent W decays are not independent. Obviously, this includes 
all Feynman diagrams in which a particle is exchanged between the production subprocess 
and one of the decay subprocesses or between the decay subprocesses. Examples for such 
manifestly non-factorizable corrections are the diagrams (a), (b), and (c) in Fig. [I]. If 
the additional exchanged particle is massive, the corresponding correction has no double 
pole for on-shell W bosons. However, if a photon is exchanged between the different sub- 
processes, this leads to a doubly-resonant contribution originating from the soft-photon 
region. This can be directly seen from the usual soft-photon approximation (SPA), which 
yields contributions proportional to the (doubly-resonant) lowest-order contribution. 

The doubly-resonant contributions can be extracted on the basis of a simple power- 
counting argument. For instance, the loop integral corresponding to diagram (c) in Fig. [1| 
is of the following form: 



d 4 g 



N(q, h 



(q 2 ~ A 2 )[(g - hf - m 2 ][(q -kJ) 2 - M%\[(q + k + Y - M^][(q + k 2 f - m 2 } 



(a) type (mf') (b) type (mf') 




Figure 1: Examples of non-factorizable photonic corrections in 0(a). The shaded blobs 
stand for all tree-level graphs contributing to e + e~ — > W + W~. Whenever Feynman dia- 
grams with intermediate would-be Goldstone bosons ± instead of W ± bosons are relevant, 
the inclusion of such graphs is implicitly understood. 



,4 N(q, kj) 

V - A 2 )(g 2 - 2qk 3 )(q 2 - 2qk_ + k 2 _ - M^)(g 2 + 2qk + + k\ - M^)(g 2 + 2qk 2 ) ' 



(2.8) 



where we have introduced an infinitesimal photon mass A to regularize the infrared (IR) 
singularity. The function N(q, fcj) involves the numerator of the Feynman integral, i.e. a 
polynomial in the momenta q and fcj, and possible further denominator factors originating 
from propagators (hidden in the blob of the diagrams) that are regular for q = and 
k± = Myj- For on-shell W bosons (k± = Myj), the integral has a quadratic IR singularity. 
For off-shell W bosons, part of the IR singularity is regularized by the off-shellness, k\ — 
M^f 7^ 0, such that the usual logarithmic IR singularity remains. Vice versa, the off-shell 
result develops a pole if either W boson becomes on shell, and is thus doubly-resonant. 
Thus, the quadratic IR singularity in the on-shell limit is characteristic of the doubly- 
resonant non-factorizable contributions. All terms that involve a factor q in the numerator 
are less IR-singular and therefore do not lead to doubly-resonant contributions and can 
be omitted. Similarly, q can be neglected in all denominator factors included in N(q, k{). 
In summary, q can be put to zero in N(q, kj) in DPA. We have checked this for various 
examples explicitly. As a consequence, we are left with only scalar integrals, and the 
non-factorizable virtual corrections are proportional to the lowest-order matrix element. 
We call the resulting approximation extended soft-photon approximation (ESPA). It differs 
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from the usual SPA only by the fact that q is not neglected in the resonant W propagators. 
In ESPA, diagram (c) in Fig. [l| gives the following contribution to the matrix element: 



M hh = ie 2 Q 2 Q 3 M Bo m 



d A q 4k 2 k 3 



(2vr) 4 (g 2 - A 2 )(g 2 - 2qk 3 )(q* + 2qk 2 ) 



x {kl-M^kl-Mlj, 

[{q-k.f-MU{q + k + ) 2 -MW 



The q 2 terms in the last four denominators are not relevant in the soft-photon limit 
and were omitted in Refs. ^, [RJ. In fact, using the above power-counting argument it 



can easily be seen that the differences of doubly-resonant contributions with and without 
these q 2 terms are non-doubly- resonant. We have chosen to keep the q 2 terms, because we 



want to use the standard techniques for the evaluation of virtual scalar integrals JTT]. In 
DPA, i.e. if we perform the limit h± — > after evaluating the integral, we should obtain 
the same result. 

In order to arrive at physical results, we have to incorporate the finite width of the 
W bosons. In DPA this can be done in at least two different ways: 

As a first possibility, we perform the integrals for zero width and afterwards put k± = 
M^f where this does not give rise to singularities. In all other places, i.e. in the resonant 
propagators and in logarithms of the form ln^ 2 , — + ie), we replace k± — M(y + ie by 
K± = kl-M^ + iM w T w- Since the width is only relevant in the on-shell limit, it is clear 
that the (physical) on-shell width has to be used. 

Alternatively, we introduce the width in the W propagators before integration. This 
has to be done with caution. If we introduce the finite width by resumming W-self-energy 
insertions, the width depends on the invariant mass of the W boson and thus on the 
integration momentum. Fortunately, the contribution we are interested in results only 
from the soft-photonic region where the virtual W bosons are almost on shell. Therefore, 
we can insert the on-shell width inside the loop integral. After performing the integral, 
we put k\ = M^r and T w = where this does not lead to singularities. In DPA this gives 
the same results as the above treatment. 

In the following we write M 2 = — iM\yPw instead of M{y in the loop integrals. It 
is always understood that M 2 and k\ are replaced by where possible after evaluation 
of the integrals. 

If we implement the width into the integrand, it is clear that only the part of the inte- 
gration region with |go| ~ Tw contributes in DPA. If | go I Tw, one of the W propagators 
must be non-resonant and the contribution becomes negligible. 

Once the width is introduced, it becomes evident that the relative error of the DPA is of 
the order of r w /scale. Let Ecm be the centre-of-mass (CM) energy and AE = E C m~ 2M w 
be the available kinetic energy of the W bosons. Then, for AE > M\y the scale is given 
by Mw, for Tw ^ AE < Mw it is given by AE and for AE < Tw it is given by Tw- This 
shows that the DPA is only sensible several Tw's above threshold. This is simply due to 
the fact that, close to threshold, the phase space where both W propagators can become 
doubly-resonant is very small, and the singly-resonant diagrams become important. 
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Figure 2: Further examples of non-factorizable photonic corrections in 0(a). 



2.3 Classification and gauge-independent definition of the non- 
factorizable doubly-resonant virtual corrections 

Manifestly non-factorizable corrections arise from photon exchange between the final states 
of the two W bosons (ff'), between initial and final state (if), and between one of the in- 
termediate resonant W bosons and the final state of the other W boson (mf ) . Examples 
for these types of corrections are shown in Fig. [I] (c), Fig. |2] (a), and Fig. [I] (a,b), respec- 
tively. In addition, there are diagrams where the photon does not couple to uniquely 
distinguishable subprocesses. These contributions can be classified into photon-exchange 
contributions between one of the intermediate resonant W bosons and the final state of 
the same W boson (mf), between the intermediate and the initial state (im), between the 
two intermediate W bosons (mm'), and within a single W-boson line, i.e. the photonic 
part of the W-self-energy corrections (mm). Diagrams contributing to these types of cor- 
rections are given in Fig. |2| (c), Fig. |^ (b), Fig. [1] (d), and Fig. ^] (d), respectively. Because 
the photon coupling to the W boson can be attributed to the decay or the production 
subprocesses, these diagrams involve both factorizable and non-factorizable corrections. 

In order to define the non-factorizable corrections, we have to specify how the factoriz- 
able contributions are split off. This should be done in such a way that the non-factorizable 
corrections become gauge- independent. In Refs. || [l(J this was reached by exploiting the 



fact that in ESPA the matrix element can be viewed as a product of the lowest-order ma- 
trix element with two conserved currents. Taking all interferences between the positively 
and the negatively charged currents arising from the outgoing W bosons and fermions 
gives a gauge-independent result. 

We have chosen a different definition of the non-factorizable corrections, which, how- 



ever, turns out to be equivalent to the one of Refs. H, 10] in DPA. Our approach has the 
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advantage of providing a clear procedure how to combine factorizable and non-factorizable 
contributions to the full 0{a) correction in DPA. Because the complete matrix element is 
gauge-independent order by order, the sum of all doubly-resonant 0(a) corrections must 
be gauge-independent. On the other hand, the factorizable doubly-resonant corrections 
can be defined by the product of gauge-independent on-shell matrix elements for W-pair 
production and W decays and the (transverse parts of the) W propagators, 



M t = £ ^{5M« + «-^-MZ7 hh MZ.7 hh (2-10) 



A +1 A_ * ^ 

•^'Born vjvi - /vl Born " /vl Born ■ /vl Born 

where 5A4 e+e ^ w+w ; SA4 W+ ^^^ 2 , and 5A4 W denote the one-loop amplitudes of the 

respective subprocesses. We can define the non-factorizable doubly-resonant corrections 
by subtracting the factorizable doubly-resonant corrections from the complete doubly- 
resonant corrections. This definition allows us to calculate the complete doubly-resonant 
corrections by simply adding the factorizable corrections, defined via the on-shell matrix 
elements, to our results. Our definition can be applied diagram by diagram. In this way, 
all diagrams that are neither manifestly factorizable nor manifestly non-factorizable can 
be split. Such diagrams receive doubly-resonant contributions from the complete range 
of the photon momentum q, and not only from the soft-photon region. This is obviously 
due to the presence of two explicit resonant propagators. However, after subtracting 
the factorizable contributions, all doubly-resonant terms that are not IR-singular in the 
on-shell limit cancel exactly, i.e. only the soft-photon region contributes. Consequently, 
also in this case q can be neglected everywhere except for the denominators that become 
IR-singular in the on-shell limit. As an example, we give the non-factorizable correction 
originating from diagram (d) of Fig. [I]: 1 



yVl nf ~ ie At Born 



d 4 g Ak+k. 



(2vr) 4 q 2 [(q + k + ) 2 - M 2 } [(q - k^) 2 - M 2 } 
d 4 q 4k + k_ 
(2tt) 4 {q 2 - \ 2 ){q 2 + 2qk + ){q 2 - 2qk_ 



:2.ir 



This example shows that the on-shell subtraction introduces additional IR singularities. If 
the IR singularities in the non-factorizable real corrections are regularized in the same way, 
they cancel in the sum. In ( p. 11 ) an infinitesimal photon mass A is used as IR regulator, 



but we have repeated the same calculation also by using a finite W-decay width as IR 
regulator instead of A, leading to the same results in the sum of virtual and real photonic 
corrections. 

We illustrate our definition of the non-factorizable corrections also for the photonic 
contribution to the W-self-energy correction [diagram (d) of Fig. 0. The non-factorizable 
part of the W + self-energy reads 



MT W+ - -ie 2 M Bmn {J 



(2vr) 4 q 2 [(q + k + ) 2 - M 2 }(k\ - M 2 ) 



1 ^Ne use the sign ~ to indicate an equality within DPA, i.e. up to non-doubly-resonant terms. 
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d 4 q 



1 



(2vr) 4 g 2 (g 2 + 2qk + ) J k 2 +=M 2 k\ - M 2 

d 4 g Ak 2 + 
(27r) 4 (g 2 - A 2 )(g 2 + 2g^ 



(2.12) 



The first integral results from the off-shell self-energy diagram, the second from the corre- 
sponding mass-renormalization term, and the third integral is the negative of the on-shell 
limit of the first two integrals. The integrals in ( p.!2j ) are UV-divergent and can be easily 
evaluated in dimensional regularization. 

The gauge independence of the non-factorizable corrections has been ensured by con- 
struction. The consistent evaluation of gauge theories requires, besides gauge independence 
of the physical matrix elements, the validity of Ward identities. It was found in Ref. f!2 
that the violation of Ward identities can lead to completely wrong predictions. The pro- 
cedure described above for extracting the non-factorizable corrections from the full matrix 
element does not lead to problems with the Ward identities that rule the gauge cancella- 
tions inside matrix elements. This is due to the fact that the non-factorizable corrections 
are proportional to the Born matrix element. Therefore, if Ward identities and gauge 
cancellations are under control in lowest order, the same is true for the non-factorizable 
corrections. 

Finally, we show how our definition of the non-factorizable corrections can be rephrased 
in terms of products of appropriately defined currents. By using 



1 



1 



(q ± k±) 2 — M 2 q 2 ± 2qk± 



1 - 



k\ 



M 2 



(q±k±) 2 -M 2 



(2.13) 



and the fact that in DPA k\ can be put to Mjy before integration in integrals that do not 
depend on M 2 , the contribution ( p. 11 ) can be expressed as 



M n{ ~ 



ie 2 M 



Born 



d 4 g 



4fc+ k. 



(2tt) 4 (q 2 - X 2 )(q 2 + 2qk + )(q 2 - 2qk. 
k 2 _ — M 2 kl- M 2 



kl-M 2 



(q + k + ) 2 - M 2 
k 2 -M 2 



(q-k_) 2 -M 2 (q + k + ) 2 - M 2 (q - k_) 2 - M 2 



(2.14) 



The other non-factorizable corrections that involve photons coupled to W bosons can 
be rewritten in a similar way. Finally, all non-factorizable virtual corrections can be cast 
into the following form: 



M 



virt 
nf 



iM 



Born 



d 4 g 



(2vr) 4 q 2 - A 2 



■ . e +e-— W+W-,u •W+-+A/2 , -c+e-^W+W" ,m -Vf-^fsU 
Jvirt Jvirt,^ ' Jvirt Jvirt,/^ 



' Jvirt Jvirt,fi 



(2.15) 



with 



■ e -r e -_> w + w - __ / 2fe +M 2/c_ M 2p_ M 2p_, 

Jvirt, u i _.9 i n_.T. ' 



Virt ^ " \q 2 + 2qk + q 2 -2qk„ q 2 - 2qp_ q 2 + 2qp^ 
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Figure 3: Example of a non-factorizable real correction. 



.w^/,/ 2 _l n 2fc lft 2k 2fl 2fc +ft \ kl-M 2 

Jvirt.u - e VI 9 , n , V2 



Vlrt '" ~~ q 2 + 2qk l ^q 2 + 2qk 2 q 2 + 2qk + J (k + + q) 2 - M 2 1 

JviTt ^ ' = " e l g3 g2 _ 2g fc 3 ~ g 2 - 2qh + q 2 ~ 2qkJ (fc_ - g) 2 - M 2 ' (2 ' 16) 



The last term in (|2.15 ) originates from the Feynman graphs shown in Fig. [I] and those 



where the final-state fermions are appropriately interchanged [interference terms (ff' 
(mf'), and (mm')]. The contributions involving the current „ _+w w contain the inter- 
ference terms (if), (mf), (im), (mm), and the remaining contributions of (mf') and (mm'). 
The contribution of the W + self-energy is given, for instance, by the product of the two 

terms involving k +fl in jvirtjT W+W ~ and Jvil^ hh - 

In DPA, the q 2 terms in the denominators of (|2.16|) can be neglected, and the cur- 
rents are conserved. The currents j^I/T^ 2 an d jXt./T^ 4 are the ones mentioned in 
Refs. H |10|| . For the virtual corrections, this shows that our definition of non-factorizable 
doubly- resonant corrections coincides with the one of Refs. || [HJ in DPA. 

2.4 Doubly-resonant real corrections 

The photonic virtual corrections discussed above are IR-singular and have to be combined 
with the corresponding real corrections in order to arrive at a sensible physical result. The 
real corrections originate from the process 

e + (p+) + e~(p_) -> W + (*4)+W~(A;'_)[ + 7(g)] 

- h(K) + fM) + f 3 (k' 3 ) + U{k' A ) + 7 (g). (2.17) 

Note that we have marked the fermion momenta k\ by primes in order to distinguish 
them from the respective momenta without real photon emission. The momenta of the 
W bosons are k' + = k[ + k' 2 and k'_ = k' 3 + k! A if the photon is emitted in the initial state 
or in the final state of the other W boson, and k + = k[ + k' 2 + q or k_ = k' 3 + k' A + q if the 
photon is emitted in the final state of the respective W boson. 

The non-factorizable corrections induced by the process ( 2.17|) arise from interferences 



between diagrams where the photon is emitted from different subprocesses. A typical 
non-factorizable contribution is shown in Fig. 0. Including the integration over the photon 
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phase space, this contribution has the following form: 




d 3 q N(q, K)5( P+ + p_ - k[ -k' 2 -k> 3 -k' 4 -q) 
2<?o [(q + k' 2 ) 2 - ml][{q + k' + ) 2 - M 2 }(k' 2 - M 2 ) 



I 



\ [(q + k' 3 ) 2 - m 2 }(k' 2 - M 2 )[(q + k'_) 2 - M 2 ] 
d 3 q N(q, fcQ£(p + + p_ - k[ -k' 2 ~k' 3 -k' 4 -q) 
2q (2qk' 2 )(2qk' + + k' 2 - M 2 ){k' 2 - M 2 ) 



9o=v /c l 2 + A2 



1 



(2.18) 



(2qk' 3 )[k' 2 - (M*) 2 ][2qkL + k' 2 - (M*) 2 ] go =V^' 



where we again use a photon mass A to regularize the IR singularities, and N(q, kQ has 
the same meaning as above. Again, the doubly-resonant contributions are characterized 
by a quadratic IR singularity for k± = M^, T w — > 0, and only soft-photon emission is 
relevant in DPA. For this reason, the W bosons are nearly on shell, and the on-shell width 
is appropriate. As for the virtual corrections, the introduction of the width before or after 
phase-space integration leads to the same results in DPA. As already indicated in ( p. 18 ), 
in the following real integrals we use M 2 with the understanding that it has to be replaced 
by My? after integration where possible. 

The aim is to integrate over the photon momentum analytically and to relate the 
fermion momenta k[ to the ones of the process without photon emission, k{. Primed and 
unprimed momenta differ by terms of the order of the photon momentum: k\ = ki + 0(qo). 
In DPA we can neglect q in N(q, /cQ, leading to the replacement N^q.k'^) — > N(0,ki). 
Moreover, we can extend the integration region for go to infinity, because large photon 
momenta yield negligible contributions in DPA. After extension of the integration region 
the integral becomes Lorentz-invariant. 

While the correction factor to the lowest-order cross section is universal in SPA for 
all observables, the correction factor is non-universal in ESPA. In order to define this 
correction factor in a unique way, one has to specify the parametrization of phase space, 
i.e. the variables that are kept fixed when the photon momentum is integrated over. This 
fact has not been addressed in the literature so far. 

Let us consider this problem in more detail. It can be traced back to the appearance 
of the photon momentum q in the 5-function for momentum conservation. In the usual 
SPA q is neglected in this 5-function, which is sensible if the exact matrix element is a 
slowly varying function of q in the vicinity of q = 0. However, in the presence of resonant 
propagators, in which q cannot be neglected, the simple omission of q in the momentum- 
conservation 5-function leads to ambiguous results: putting q = in the ^-function and 
identifying k[ with ki in ( |2.18| ) yields 




1 



(2.19) 



(2qk 3 )[k 2 + - (M*) 2 ][2qk_ + hi - 
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On the other hand, eliminating k' + in the denominator of (|2.18 ) with the help of the 8 



function, putting q = in the ^-function, restoring k' + with the modified 5-function, and 
setting k[ — > hi results in 



1 



d 3 q N(0, kj)6(p + + p- - fei - k 2 - k 3 - fc 4 ) 
2q (2qk 2 ) {k\ - M 2 ){kl - M 2 ) 

(2.20) 



(2qk 3 )[-2qk + + k\ - (M*) 2 ][2qk„ + kl - (M*) 



<jty 



Both expressions differ by a doubly-resonant contribution. The difference is in general 
confined to the W propagators and originates from the fact that not only soft photons but 
also photons with energies of the order of \k± — M^|/M w or, after the inclusion of the 
finite width, of order Tw contribute in DPA. Since photons with finite energies contribute, 
it is evident that the integral over the photon momentum depends on the choice of the 
phase-space variables that are kept fixed. 

As a consequence, one has to choose a definite parametrization of phase space and to 
exploit the 5-function carefully, in order to define the non-factorizable corrections uniquely. 
For instance, if the vector k' + = k[ + k' 2 is kept fixed, the alternative ( 2.20j ) is excluded. 



However, because of momentum conservation, not all external momenta can be kept fixed 
independently. 

It is, however, possible to keep, for instance, the invariant masses of the final-state 
fermion pairs k' 2 = (k[ + k' 2 ) 2 and k' 2 = (k' 3 + k' 4 ) 2 fixed when integrating over the photon 
momentum. If we require {k[ + k 2 ) 2 = k' 2 = k 2 ^ = + A^) 2 , we obtain for the denominator 
of the W + boson 

(q + k '+f -M 2 = 2qk' + + k'l-M 2 = 2qk + + k\ - M 2 + C(g 2 ) 

= (q + k + ) 2 -M 2 + 0(q 2 ), (2.21) 

where k\ = hi + O(q ) was used. Based on the power-counting argument given above, the 
terms of order q 2 can be neglected in DPA, and we find 

(q + k' + ) 2 -M 2 ~ (q + k + ) 2 -M 2 . (2.22) 

If we choose to eliminate k' + , as done in the derivation of (2.20), we find, on the other 
hand, 

(q + k' + ) 2 — M 2 = (p + + p_ - k'_) 2 — M 2 = (p + + p_ - k_f — M 2 + O(q ) 

= k\- M 2 + O(q ). (2.23) 



The O(qo) terms are relevant in DPA [and in fact given by (|2.21| )]. As a consequence, ( |2.20|) 



is not correct if we choose to fix k' 2 = k\ when integrating over the photon momentum. 
For fixed k' 2 = k\, (|2.21|) leads to the unique result (|2.19| ) for the W + propagator in DPA, 



independently of the other phase-space parameters. If we choose, on the other hand, to fix 
Arj_ = Arj_, which corresponds to a different definition of the invariant mass of the W + boson, 
we obtain 

(q + k' + ) 2 — M 2 = k 2 + — M 2 , (2.24) 
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and thus ( |2.20|) instead of (|2.19|) . Consequently, the different approaches ( |2.19| ) and (|2.20|) 
correspond to different definitions of the invariant mass of the W + boson which decays into 
the fermion pair (fi, f 2 ). In order to define the DPA for real radiation, one has to specify 
at least the definition of the invariant masses of the W bosons that are kept fixed. In the 
following we always fix k\ = (k' l +k' 2 ) 2 = (&1+&2) 2 and k 2 _ = (k^+k'^) 2 = (ks+k^) 2 , as it was 
also implicitly done in Refs. || |9|. [l(J . Once the invariant masses of the W bosons are fixed 
in this way, the resulting formulae for the non-factorizable corrections hold independently 
of the choice of all other phase-space variables. 

We stress that the results obtained within this parametrization of phase space differ 
from those in other parametrizations by doubly-resonant corrections. As already indicated 
in the introduction, in an experimentally more realistic approach the invariant masses of 
W bosons are identified with invariant masses of jet pairs, which also include part of the 
photon radiation. Since this situation can only be described with Monte Carlo programs, 
our results (as well as those of Refs. p, || |10j) should be regarded as an estimate of the 
non-factorizable corrections. Moreover, our analytic results provide benchmarks for such 
future Monte Carlo calculations. 



2.5 Classification and gauge-independent definition of the non- 
factorizable doubly-resonant real corrections 

The doubly-resonant real corrections can be classified in exactly the same way as the 
virtual corrections. For each virtual diagram there is exactly a real contribution, which 
we denote in the same way, e.g. ff' refers to all interferences where the photon is emitted 
by two fermions corresponding to the two different W bosons. 

As in the case of the virtual corrections, one has to define the non-factorizable real 
corrections in a gauge-independent way. To this end, we proceed analogously and define 
the non-factorizable real corrections as the difference of the complete real corrections and 
the factorizable real corrections in DPA. The factorizable corrections are defined by the 
products of the matrix elements for on-shell W-pair production and decay with additional 
photon emission. There are three contributions to the factorizable real corrections, one 
where the photon is radiated during the production process and two where it is emitted 
during one of the W-boson decays. The corresponding matrix elements read 



A! real, 1 = 



A4 re aJ,2 = 



A+,A_ 



e+e-^W+W-7 A/f W+^/i/ 2 Ayf W "^/3/4 
^'Bom yvt Born ^'Born 

[k\ - M 2 )(k 2 - M 2 ) 

A^e+e-^W+W" A / /W+-»/i/ 2 7 \A^'^hh 
^'Bom -^Born -^Born 

[{k + + q) 2 — M 2 ](k 2 — M 2 ) 



A A V* Born J l Born J l Born (r, r>r\ 

-WreaU = 2^ 772 T\J2\\(1 1 ^2 a721 ' \i-iO) 



(k\ — M 2 )[(k^ + q) 2 — M 2 ] 



in analogy to ( |2.7| ). Note that the matrix elements A"l r eai,2 and A^reai^ involve an explicit 
g-dependent propagator. The factorizable corrections are given by the squares of these 
three matrix elements, and include by definition no interferences between them. 
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Figure 4: Real bremsstrahlung diagram containing non-factorizable and factorizable con- 
tributions. 



As an example for the extraction of non-factorizable corrections from real diagrams that 
involve both factorizable and non-factorizable corrections, we consider the contribution of 
the diagram in Fig. After subtraction of the factorizable contribution, which originates 
from | A1 r eai,i 1 2 > it gives rise to the following correction factor to the square lA^Boml 2 of the 
lowest-order matrix element fl2.7p : 

d 3 q 



W+W" 
real 



(27r)32g 



2 Re 



4k+k_ 



[(k+ + q) 2 - M 2 ][(A;_ + q) 2 - (M*) 2 ] 







2qk + 2qk_ _ 





<hv- 



Vq 2 +A 2 



(2.26) 



/ \2 



Note that the form of the correction factor is only correct for fixed {k^ + k^) 2 and (k' 3 + k 
For other conventions the off-shell contribution changes, whereas the on-shell contribution 
stays the same. Using the relations (|2.13 ) for q 2 = 0, in DPA we can rewrite Q2.26Q as 



rW+W" 
real 



d 3 



q 



(2n) 3 2q 
kl- 



2 Re 
- M 2 



(qk+)(qk-) 



k 2 



M 2 



+ 



k 2 _ - (M*) 2 



(A;+ + g) 2 -M 2 (Jfe_ + q) 



*\2 



k_ - (M*) 



(k+ + q) 2 - M 2 (jfc_ + q) 2 - (M*) 2 



(M*) 2 

(2.27) 



In the same way, all other contributions that originate from a photon coupled to a 
W boson can be rewritten such that the complete real non-factorizable corrections can 
finally be expressed as the following correction factor to the lowest-order cross section, 

d3( l ot^_ r,e+e-- > W+W-,M/ •W+-*/i/ax* , -e+e-^W+W-^/ |W--»/ 3 / 4 \* 

\Jrea\,fj, ) ' Jreal wreal,/i / 



5, 



rcal,nf 



(2vr) 3 2g 



2 Re 



•e'e 
Jreal 



with the currents 



■e+e-^W+W" 
Jveal,^, 



J real,^t 



■W-^/3/4 
Jreal,/i 



1 ■ vv 1 - 
*?real 



\qk+ 



\3 



qk 



real,^t 



\L _|_ 



n fcl/x n ^2 M 

vn — V2— r- 



qh 



Q^E _| 

qk± qk 



qk + J (k. 
k_ 



M 2 



k 2 — M 2 



{k_ + q) 2 -M 2 ' 



(2.28) 



(2.29) 
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The factor ( p. 28 ) for the non-factorizable correction can be viewed as the interference 



contributions in the square of the matrix element (e denotes the polarization vector of the 
photon) 

.^l real — -^l Born t {JreaX^ + Jreal,/i + Jreal,^ J' (^.<3UJ 

which is just the sum of the three matrix elements Aureal,,, i = 1,2,3, in ESPA (includ- 
ing radiation from the external fermions and the internal W bosons). The respective 
squares of these three contributions correspond to the factorizable corrections. Note that 
-^Bom^ W+W £M Jreai^ W+W i s the soft-photon matrix element for on-shell W-pair pro- 
duction. Similarly, M-Born £fl JreaijT'^ 2 an d ^B^rrT^^^^Jreai,^^ 3 ^ 4 correspond to the 
soft-photon matrix elements for the decays of the on-shell W bosons, apart from the extra 
factors (k± — M 2 )/[(k± + q) 2 — M 2 ]. These factors result from the definition of the lowest- 
order matrix element in terms of k± and from the fact that we impose k± = k± and not 
12 — p 

ru_|_ — ru_|_. 

Obviously, the currents ( |2.29|) are conserved, and the corresponding Ward identity for 
the U(l) em symmetry of the emitted photon is fulfilled. 

3 Calculation of scalar integrals 

In this section we set our conventions for the scalar integrals. We describe the reduction 
of the virtual and real 5-point functions to 4-point functions and indicate how the scalar 
integrals were evaluated. More details and the explicit results for the scalar integrals can 
be found in the appendix. 

3.1 Reduction of 5-point functions 
3.1.1 Reduction of virtual 5-point functions 

The virtual 3-, 4-, and 5-point functions are defined as 

C (pi,P2,m ,m 1 ,m 2 ) = 7^ / d 4 q 1 , 
D{o^}{pi,P2,P3,m ,m 1 ,m 2 ,m 3 ) = — d q , 

17T J 1\ q±\ i±\ 2-LV 3 

E{o tfl }ipi,P2,P3,P4,mQ,mi,m2,m 3 ,m 4 ) = f d A q ^l^l AT , (3-1) 

with the denominator factors 

N = q 2 - m 2 + ie, Ni = {q + Pif - m 2 + ie, % = 1, . . . , 4, (3.2) 
where ie (e > 0) denotes an infinitesimal imaginary part. 
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The reduction of the virtual 5-point function to 4-point functions is based on the fact 
that in four dimensions the integration momentum depends linearly on the four external 
momenta Pi ||13||. This gives rise to the identity 



2q 2 2qp\ . . . 2qp 4 
2p ± q 2p\ ... 2p x p 4 



27V + Foo 2qpi ••• 2gp 4 
Ni - N + Y 10 - Y 00 2p\ ... 2 Pl p 4 



2p 4 q 2p 4 p x ... 2p\ N 4 -N + Y 40 - Y 00 2p 4Pl ... 2p\ 



(3.3) 



with 



Y 00 = 2m 2 , Y i0 = Y 0i = m 2 +m 2 -p 2 , F„ = m 2 +m 2 -(pi-pj) 2 , i, j = 1, 2, 3, 4. (3.4) 
Dividing this equation by iVoiVi ■ ■ ■ N 4 and integrating over d 4 g yields 



= i(dV 1 



ivr 2 7 a iV yYi • • • N 4 



2N + Y 00 2q Pl ... 2qp 4 

Ni - Nq + Yw -Y 00 2p\ ... 2 Pl p 4 

N 4 - N + F 40 - ^oo 2p 4 Pi • • • 2p 2 



(3.5) 



Expanding the determinant along the first column, we obtain 



2D (0)+Y 00 E 



2p x p\ . . . 2p x p 4 



2p 4 p 1 . . . 2p 4 p 4 

+ J2(- 1 ) k { D ^ k )-[ D ^)+P^ D o(°)} +P^D (0) + (Y ko -Y 00 )E fl } 



k=l 



X 



2p1 
2pipi 



2p£ 
2pip 4 



2p k -iPi ■ ■ ■ 2p k _ x p 4 
2pk+iPi ■ ■ ■ tyk+m 

2p 4 p x . . . 2p 4 p 4 



where D (k) denotes the 4-point function that is obtained from the 5-point function E 
by omitting the kth propagator Nj^ . The terms involving p 4fJ ,D (0) have been added for 
later convenience. 



(3.6) 
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All integrals in (3^6) are UV-finite and Lorentz-covariant. Therefore, the vector inte- 
grals possess the following decompositions 



A,(0)+jvD o (0) 



E iPi» 

8=1 



fc = l,2,3,4, 



i=l 



EA(0)(pi-p 4 )/ 



(3.7) 



The last decomposition becomes obvious after performing a shift q — > g— p 4 in the integral. 
From (|3.7| ) it follows immediately that the terms in ( |3.6| ) that involve D^(k) drop out when 
multiplied with the determinants, because the resulting determinants vanish. Similarly, 
D^(0) +p4^D (0) vanishes after summation over k. Finally, the termp 4 ^D (0) contributes 
only for k = 4, where it can be combined with the first term in ( |3.6|) . Rewriting the 
resulting equation as a determinant and reinserting the explicit form of the tensor integrals 
leads to 

iVo + Yqo 2gpi . . . 2qp A 



Na 



Yio - Yqo 2pxpx . . . 2p 1 p 4 



F 40 - Y 00 2p 4 Pi • • • 2p 4 p 4 



(3.8) 



Using 



2pipj = Yij - Y i0 
2q Pj = Nj - N 



Y 



0.7 



Y 



Yoo, 
Yqo, 



(3.9) 



adding the first column to each of the other columns, and then enlarging the determinant 
by one column and one row, this can be written as 



1 F)0 

D (0) + Y 00 E 
iio — Yqq 



Yqa 

D (A) + Y 04 E 

Yu — Yq4 



Yaq — Y( 



00 



Yaa — Yr 



01 



(3.10) 
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Equation (|3.10|) is equivalent to 



-E D (0) D„(l) D (2) D„(3) D„(4) 



1 


Yoo 


Y01 


Y02 


Y03 


Y04 


1 


Y w 


Y u 


Y 12 


Y 13 


Y u 


1 


Y20 


Y 21 


Y22 


Y23 


Y-2A 


1 


Y30 


Y 3 i 


Y 32 


Y33 


>34 


1 


Y40 


Kn 


Y42 


Y43 


I44 



(3.11) 



which expresses the scalar 5-point function Eq in terms of five scalar 4-point functions 



(3.12) 



where Y = (Yij), and Yi is obtained from Y by replacing all entries in the ith column 
with 1. 

In the special case of an infrared singular 5-point function we have 



Yqo — 2A — >• 0, Yq\ 
and the determinants fulfil the relations 



m i - pi = °) 



Y 



04 



m\-p\ = 0, 



(3.13) 



= det(y) - Y 02 det(y 2 ) - Y 03 det(Y 3 ), (3.14) 

= (y 24 - y 02 ) det(y) + y 02 Yi 4 det(yo + (y 02 y 3 4 - y 03 y 2 4) det(y 3 ) + y 02 y 4 4 det(y 4 ), 
= (y 31 - y 03 ) det(y) + y 03 y n det(y) + (y 03 y 21 - y 02 y 31 ) det(y 2 ) + y 03 y 4 i det(y 4 ), 



which allow the simplification of ( 3.12 ). 



3.1.2 Reduction of real 5-point functions 

The real 3-, 4-, and 5-point functions are defined as 

C Q (pi,p2,\,m 1 ,m 2 ) 
Do(Pi>P2,P3, X,mi,m 2 ,m 3 ) 
£o{Pi,P2,Pb, Pi, A, m 1 ,m 2 ,m 3 , m 4 ) 



d 3 q 1 



7T Jq <AE 2q N{N^ 

2 fd 3 q 1 



vr J 2q N[N^ 
2 rd 3 q 1 



) 

9o=vV+ A2 



with 



tt7 2g NIN2N3N4 



i = l,...,4. 



(3.15) 



(3.16) 



The shift of the integration boundary of go to infinity leads to an artificial UV divergence 
in the 3-point function Co, which is regularized by an energy cutoff AE —>■ 00. In the 
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Re{g } 



Rc{9 } = |q|-e 



Figure 5: Illustration of the contour C of fl3.17| ) in the complex go plane. The open circles 
indicate the "particle poles" located at go(pi) = (2qp« — p 2 + mf) j (2p i0 ) . 



following we only need differences of 3-point functions that are independent of AE and 
Lorentz- invar iant . 

Because of the appearance of UV-singular integrals in intermediate steps, the reasoning 
of the previous section cannot directly be applied to Sq. Therefore, we rewrite the real 
5-point function as an integral over a closed anticlockwise contour C in the go plane and 
introduce a Lorentz-invariant UV regulator A: 



£o(Pi,P2,P3,P4, X,m 1 ,m 2 ,m 3 ,m^ = lim — 

A^oo \7T Z JC 



d 4 g 



1 



-A 2 



N&Ni ■ ■ ■ Ni g 2 - A 2 



(3.17) 



with 



jV = g 2 _A 2 . (3.18) 

The contour C is chosen such that it includes the poles at g = y/q 2 + X 2 and go = 
Vq 2 + A 2 , but none else (see Fig. ||)[]. 

The integral ( p. 17 ) can be reduced similarly to the virtual 5-point function. Owing to 
the different propagators (|3.3|) leads to 



If 1 
= lim — / d 4 g --7—^— 

A^oo ITT 2 JC NqN[ ■ 



-A 2 



Ni q 2 - A 2 



2N' Q 2q Pl 
N[ + Y 10 2p 2 

N' A + F 40 2p4Pi 



2gp 4 
2pip 4 



2 P l 



(3.19) 



2 It is straightforward to check that the naive power counting for the UV behaviour in ( 3.17 ) is valid 
for time-like momenta pi ; light-like pi can be treated as a limiting case. The basic idea for the proof is 
to deform the contour C to the vertical line Re{go} = | q | — e with a small e > 0, which is allowed for 
sufficiently large |q|, more precisely when all particle poles appear left from the line Re{go} = |q| — e - 
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instead of ( |3.5|) , with Y^ from ( |3.4|) , and A 2 can be set to zero in all Yy, in particular we 
have Y 0Q = 2 A 2 — > 0. After expanding the determinant along the first column and using 
the Lorentz decompositions of the integrals where N' k is cancelled, we see that these terms 
vanish, and we are left with 



A lim 



d 4 g 



-A 2 



N' Q N[ ■■■N' A q 2 -h? 



2N 2q Pl 
Y w 2p\ 



Y40 2p 4 p 1 



2 P l 



(3.20) 



Using 



%PiPj — Yij — Yio — Yoj, 
2qp j = N' j + Y 0j , 



(3.21) 



adding the first column to the other columns and extending the determinant leads to 



= lim 



A— >oo 17T Z 



d 4 g 



-A 2 



N^N{ ■ ■ ■ Niq 2 - A 2 



1 Y 01 

2N N[ + 2N' Q + F 01 

Y\q Yn — Yq\ 



Yq4 

N' 4 + 2N + F 04 
Yu — ^04 



Ya 



10 



yki - Y 



01 



Y44 — Yr 



04 



(3.22) 

Subtracting the first row from the second, adding the first row to the other rows, and 
exchanging the first two rows, we arrive at 



= lim 



d 4 g 



-A 2 



A-00 ivr 2 Jc * N^N[ ■ ■ ■ Ni q 2 - A 2 



-1 2A^ N[ + 2N' 



1 

1 ^10 



Y n 



1 Y 



10 



Ya 



41 



K + 

YqA 

Y 14 



Ya 



44 



(3.23) 



Now we perform the contour integral using power counting for A — > 00 (see footnote ||). In 
the contribution of the pole at go — Vq 2 + ^ 2 we have Nq = in the numerator; the limit 
A — > 00 can be trivially taken and reproduces usual bremsstrahlung integrals, as defined 
in ( |3.15| ). In the contribution of the pole at go — V^L 2 + A 2 the term containing Nq in the 
numerator survives and will be calculated below, but all other terms vanish after taking 
the limit A — > 00. 
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Thus, we find 

-£ V Q (0) V (1) + V q (0) V (2)+V (0) V (3)+V (0) V (A)+V (0) 



1 


Yoo 


Y01 


Y02 


Y03 


Y 4 


1 


Y 10 


Y n 


Y 12 


Y 13 


Y u 


1 


Y20 


Y 2 i 


Y22 


Y23 


Y 2 4 


1 


Y30 


Y 31 


Y32 


Y33 


Y34 


1 


Y40 


Y 41 


Y42 


Y43 


Y44 



or 



Here 



Sn 



det(y) 
V (0) 



det(Fo) £o(0) + £ det(Yi) \v (i) + V (Q) 



i=i 



lim 



d A q 



-A 2 



(3.24) 
(3.25) 

(3.26) 



A-ooivr 2 ./c * JVj---JV£g a -A 2 ' 
and the 4-point bremsstrahlung integrals T>o(i), i = 1,2,3,4, result from Sq by omitting 
the ith denominator N[. The result (|3.24 ) differs from ( 3.10|) only by the extra X> (0)'s 
added to the £> W S - 

The integral T>q(0) stems from the terms involving Nq in the numerator in (|3.23| ) and 
can be expressed as follows: 



^o(O) 



in 2 Jc 



d A q' 



-1 



2q'pi ■ ■ ■ 2q'p i q 



12 



(3.27) 



where the contour C surrounds q' 
yields 



q' + 1. Performing the contour integral over dq' 



f> (0) 



d 3 q' 1 
2q' Q 2q'pt ■ ■ ■ 2q'p 4 



q' 2 +l 



(3.28) 



Now the vector q' is time-like. Since also the vectors Pi are time-like (or at least light-like), 
the scalar products q'pi cannot become zero. After redefining the momenta, 



Pi = OiVi 



±1. 



Pio > 0, 



(3.29) 



and extracting the signs at, this integral can be evaluated by a Feynman-parameter rep- 
resentation and momentum integration in polar coordinates resulting in: 



£o(0) 



f°° / ^ 

-Gia 2 <J3<J4 / dxi dx 2 dx3 dx 4 5 ( 1 — V] 
Jo - \ \ j 



i=l 



(3.30) 



This is just the Feynman-parameter representation of a virtual 4-point function such that 
we finally obtain 



X> Q (0) = -aia 2 (T3(T A DQ (p 2 - p u p 3 - p u p 4 - Pi, Jpj, Jpj, \JPh VPl 



(3.31) 
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3.1.3 Explicit reduction of the virtual 5-point function for the photon ex- 
change between / 2 and / 3 

For the photon exchange between f 2 and / 3 the following scalar integrals are relevant: 



Eq 


= E ( 


-fa, 


-fa,k + , fa, A, m 3 , M, M, m 2 ] 


D (0) 


= D ( 


-fa, 


k + + fa, fa + fa, 0, M, M, 0), 


Do{l) 


= D ( 


-fa 


,k+,fa,0,M, M,m 2 ), 


D (2) 


= D { 


-h, 


k + ,k 2 ,\,m 3 ,M, m 2 ), 


A) (3) 


= D ( 


-fa, 


-fa,k 2 ,X,m 3 ,M, m 2 ), 


A) (4) 


= A( 


-fa, 


-fa,k + ,0,m 3 ,M, M). 



(3.32) 



Since we neglect the external fermion masses, the last two relations ( |3.14|) simplify to 








K + (s 23 + s 24 ) + K_Ml det(Y 3 



[s 23 + s 24 ) det(y) + K_s 23 det(Fi) 

[s 13 + s 23 ) det(Y) + K +S23 det(Y 4 ) - [K-(s 13 + s 23 ) + K + M 2 \ det(Y 2 ). (3.33) 
These relations allow us to eliminate det(Yi) and det(Y 4 ) from (|3.12|) , resulting in: 

det(Yo) 



E (-fa, -fa, k + , fa, X, m 3 , M, M, m 2 
det(y 3 ) 



det(Y) 



A(o) 



+ 



det (Y)K_s 23 



+ 



{[K+(s 23 + s 24 ) + K-Ml]D (l) + iT_ S23 A(3)} 
^^—{lK_(s 13 + s 23 ) + K + M 2 ]D (4) + K + s 23 D (2)} 



det(Y)K + s 23 



S13 + S23 



A(4) 



if+S 23 

The matrix Y reads 

^ -K_ 

* M 2 

Y 



S23 + S24 
Ks 23 



DM). 



(3.34) 



'-K, 



Sl3 - S23 



* * 2M 2 



2M 2 



* 
* 



2M 2 
* 





-S23 

— s 23 — s 2i) 

M 2 




(3.35) 



/ 



Neglecting terms that do not contribute to the correction factor in DPA, the corresponding 
determinants are given by 

det(y) ~ 2s 23 \K + K„s u s 23 - (K+M 2 , + K.sy^K^M 2 , + K + s 24 ) 



det(y ) 
det(y) 



-K 



W! 



^4 



^w( S 23 — s 24) + (^23 + S 2 4)(Si3S 2 4 — ^14%) 

w(-< 



K-M™{—ML + 2si 3 s 23 + si 3 s 24 + S14S23 
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det(F 2 



-s 23 



K + (M^ + s 13 s 24 - s 14 s 23 ) + 2K_M^s 



13 



det(F 3 ) = det(r 2 )| x+ ^_ iSi3 _ 24 , 
det(y 4 ) = det(r 1 )lx + ^_,s 13 ^ 24 > 
where the shorthand % is defined in Appendix ^. 



(3.36) 



3.1.4 Explicit reduction of the real 5-point function for the photon exchange 
between f 2 and f 3 

In this case the integrals appearing in (|3.25|) read 



£ 


= £ (h, fc_ 


, k+, k 2 , A, m 3 , M*, M, m 2 ), 


f>o(0) 


= -D (k 4 , 


k + -k 3 ,k 2 -k 3 ,0 } M_,M + ,0) } 


^o(l) 


= V (k_,k 


+,k 2 ,0,M*,M,m 2 ), 


V (2) 


= V (k 3 ,k- 


h ,k 2 ,X,m 3 ,M, m 2 ), 


V (3) 


= T> (k 3 , k. 


-,k 2 , \,m 3 ,M*,m 2 ), 




= V (k 3 ,k. 


_,fc + ,0,m 3 ,M*,M). 



(3.37) 

In analogy to the virtual 5-point function, we can express the real 5-point function with 
the help of the relations (we denote the matrix Y for the real 5-point function with a 
prime, in order to distinguish it from the one for the virtual 5-point function): 

= (s 23 + s 24 ) det(y') + K*_s 23 det(YD - [K + (s 23 + s 24 ) - ICM 2 + ] det(r 3 / ), 

= (s 13 + s 23 ) det(F') + K + s 23 det(F 4 ') - \k*_(s 13 + s 23 ) - K + M 2 } det(Y 2 '), (3.38) 



by 

£ (k 3 , k-, k + , k 2 , A, m 3 , AT, M, m 2 
det(Y 3 ') 



det( l^o(0) 



+ 



det(Y')Kls 23 

det(F 2 Q 
det(Y')K + s 23 

$13 + s 23 



det(y 

K + (s 23 + s 24 ) - KIM 2 ] fo (l) + Po(0)l + K*ls 23 \v (3) + P o (0) 



K*_(s 13 + s 23 )-K+M_ 

$23 + S 24 - 



Po(4)+P (0) 



K + s 23 

For the matrix Y' we find 

^ -K*_ 
* (M*) 2 
Y' - * * 2(M*) 2 (-2K 

1 * * * 



Kls 23 



P o (l)+Po(0) 



K +S23 [P (2)+D (0) 

(3.39) 



-K+ + s 13 + s 23 ) 
2K*_ + s- (M*) 2 
2M 2 
* 



M 2 





s 23 

+ s 23 + s 24 
M 2 




\ 



(3.40) 
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Replacing —K*_ by K_ and (M*) 2 by M 2 and multiplying the second and third columns 
and rows by —1, this becomes equal to (|3.35| ) in DPA. 

In DPA, Z> (0) can be neglected in the terms V (i) + Z> (0) in ( pT25|) and ( ggg ), and 



the reduction of the real 5-point function becomes algebraically identical to the reduction 
of the virtual 5-point function, apart from the differences in signs of some momenta. As 
a consequence, the results for the virtual corrections can be translated to the real case if 
we substitute K_ — > — K*_ in all algebraic factors such as the determinants and E — > £ , 

D (0) -> P (0), £> (1) -> -^o(l), A)(2) -> -X> (2), A)(3) -> £»o(3) and D (4) -> X> (4). 
In particular, the determinants are related by 



det(Y') ~ +det(F) 

det(y i / ) det(yi) 

det(y 3 ') ~ +det(y 3 ; 



K- 



-K* 



K-^-K* 



det(y ') ~+det(y ) 

det(y 2 ') ~ -det(y 2 ; 
det(y 4 ') ~ +det(y 4 ; 



~ —K 



K- 



K- 



w ■ 



-K* 



(3.41; 



3.2 Calculation of 4- and 3-point functions 



The scalar loop integrals have been evaluated following the methods of Ref. |H |. Our 



explicit results are listed in Appendix |C.lj . For vanishing W-boson width they agree with 
the general results of Ref. [I1J . For finite W-boson width the virtual 4-point functions are 
in agreement with those of Refs. 0, 10] in DPA. This shows explicitly that the q 2 terms 
in the W-boson and fermion propagators are irrelevant in DPA. 

An evaluation of the bremsstrahlung integrals, which follows closely the techniques for 
calculating loop integrals, is sketched in Appendix [FJ. The final results in DPA are listed 
in Appendix p.2j , and the 4-point functions agree with those of Refs. || |10 1 . We have 
analytically reproduced all exact results for the occurring bremsstrahlung 4- and 3-point 
integrals by independent methods. In addition, we have evaluated the IR-finite integrals 
T> (1), X ) o(4), and £ — V (3)/K + by a direct Monte Carlo integration over the photon 
momentum, yielding perfect agreement with our exact analytical results for these integrals. 
Note that this, in particular, checks the reduction of the bremsstrahlung 5-point function 
described in the previous section. 

Our results for the 3-point functions cannot directly be compared with those of 
Refs. 0, [TIJ, because different approaches have been used. While our results are IR- 
singular owing to the subtracted on-shell integrals, the results of Refs. || [10[ are artifi- 
cially UV-singular owing to the neglect of q 2 in the W propagators. However, when adding 
the real and virtual 3-point functions the two results agree. This confirms that our defi- 
nition of the non-factorizable corrections is equivalent to the one of Refs. [§, in DPA. 
Thus, it turns out that in DPA the subtraction of the on-shell contribution is effectively 
equivalent to the neglect of the q 2 terms in all but the photon propagators. 



4 Analytic results for the non-factorizable corrections 
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4.1 General properties of non-factorizable corrections 

In Ref. || it was shown from the integral representation that the non-factorizable correc- 
tions associated with photon exchange between initial and final state vanish in DPA. This 
was confirmed in Refs. || [K|. Via explicit evaluation of all integrals we have checked that 
the cancellation between virtual and real integrals takes place diagram by diagram once 
the factorizable contributions are subtracted. In this way all interference terms (if), (mf), 
(im), and (mm) drop out. Examples for the virtual Feynman diagrams contributing to 
these types of corrections are shown in Fig. |2|. 

The only non- vanishing non-factorizable corrections are due to the contributions (ff'), 
(mf), and (mm'). The corresponding virtual diagrams are shown in Fig. |], apart from 
permutations of the final-state fermions. Two of the corresponding real diagrams are 
pictured in Figs. |3| and f|. Since these corrections depend only on s-channel invariants, the 
non-factorizable corrections are independent of the production angle of the W bosons, as 
was also pointed out in Refs. || [10| . 



4.2 Generic form of the correction factor 

The non-factorizable corrections da n f to the fully differential lowest-order cross-section 
duBorn resulting from the matrix element ( |2.7| ) take the form of a correction factor to the 
lowest-order cross-section: 

ddnf = 5 n f dffBorn- (4-1) 

Upon splitting the contributions that result from photons coupled to the W bosons ac- 
cording to 1 = Qw+ = Qi — Q2 and 1 = — Qw- = — Q3 into contributions associated 
with definite final-state fermions, the complete correction factor to the lowest-order cross- 
section can be written as 



<w = E E (-1) 

o=l,2 b=3,4 



-6+1 



a 



QaQb - Re{A(k + , k a ; k b )}. 
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(4.2) 



In the following, only A = A(k + , k_, k%) is given; the other terms follow by obvious 
substitutions. As discussed above, A gets contributions from intermediate-intermediate 
(A mm /), intermediate-final (A mf /), and final-final (A ff /) interactions: 



A = A mm , + A 



mf 



(4.3) 



The quantity (A mm / is independent of the final-state fermions. The individual contribu- 
tions read 



A, 



(2M W - s)\ C (k + , -fc_, 0, M, M) - C (k + , -k_,X, M w , M, 



C (k+, k_, 0, M, M*) + C {k + , k-, A, M w , M, 



A m f ~ -(s 23 + s 24 ) K + D (l) - K + V Q (1] 



vSl3 + S23) 



A 



ff 



-K,s 



2:; 



K_E - K*l£ 



(4.4) 

(4.5) 
(4.6) 
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with the arguments of the 5- and 4-point functions as defined in ( |3.32| ) and ( |3.37|) . 

The sum A mf ' + Aff can be simplified by inserting the decompositions of the 5-point 
functions (ggp and (|]3|). In DPA this leads to 



A mf / + Aff' 



+ 



det(F) w det(Y') 

K+ d ^ 3) {[*"+(**. + s 2A ) + K„Ml\D Q {l) + K.s 23 D (3)} 

^ de d t 6 ( y[ 2) {[K-(s 13 + s 23 ) + K + M^}D (A) + K + s 23 D (2)} 
K ^Y^) ] { [K+(S23 + S24) " ^- M w]^o(l) + K*_s 23 V (3)} 



+ { [K1{Sl3 + S23) " K + M ^ Vo ^ + K + s ^o(2)} . (4.7) 

Note that A mf / is exactly cancelled by the contributions of the last two terms in ( |3.34| ) 



and flQ9l) . 

Inserting the expressions for the scalar integrals from Appendix |C| and using the first 
of the relations ( |3.14|) , various terms, notably all IR divergences and mass-singular loga- 
rithms, cancel between the real and virtual corrections, and in DPA we are left with 



A mm ' ~ 



s(3- 



w 



Li 2 \ -g-i&N J + £i2[j^-,x w j + Ci 2 \- — ,x w 



- Ci 2 \ -— ,i w + 27riln' 



K 4 



+ imaginary parts, (4.8) 



A m f + A 



ff 



K + K_s 23 det(Y ) ^ + det(Y 3 ) F K_det(Y 2 ) 

-1J0[V) ; r 3 r7T7 ^ 1 2 



det(y) 



det(Y) 



K + K*_s 23 det(Y>) ~ KtjetQg) 



det(y) 
+ imaginary parts, 



det(y) 



det(y) 
K*_ det(y 2 / ) 
det(y) 



^2 



(4.9) 



with .Do(O) and Po(0) given in ( |C3| ) and ( |C.10| ), respectively, and 



7-1 ~. (K+ S23 + s 2i . 



1 w 



E 

T=±l 



VA- / V S23 + S24 



+ ie, x 



Li 2 ( +2\n(-^l- - ie) hi(~ 

s 23 ) \ Ml, ) V ^ 2 



ir 



, 2 / S 23 + s 24 



T 3 



K-^,-K* 



+ 2vri 



2 l n( l_^t^ + ^_ ln / 1+ A' 



K*Lx w 



In 1 + 



■§23 + S 24 



In 



iK* 



S23 + s 24 



•ra — F 3 
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To 



K- 



-K* 



+ 27ri<Mn 1 + 



lnfl + XwM w 



Sl3 + S23 



+ In 



iM4(si 3 + S23) 



(4.10) 



The variables ftw, xw and the dilogarithms Li 2 , £^2 are defined in Appendix |A|. 

The above results contain logarithms and dilogarithms the arguments of which depend 
on K±. It is interesting to see whether enhanced logarithms of the form ln(K±/M^j) 
appear in the limits K± — * 0. It turns out that such logarithms are absent from the 
non-factorizable corrections, irrespective of the ratio in which the two limits K± — > are 
realized. 



4.3 Inclusion of the exact off-shell Coulomb singularity 

For non-relativistic W bosons, i.e. for a small on-shell velocity /3\y, the long range of the 
Coulomb interaction leads to a large radiative correction, known as the Coulomb singular- 
ity. For on-shell W bosons, this correction behaves like l//3\y near threshold, but including 
the instability of the W bosons the long-range interaction is effectively truncated, and the 
l//?w singularity is regularized. Therefore, for realistic predictions in the threshold re- 
gion, the on-shell Coulomb singularity should be replaced by the corresponding off-shell 
correction. The precise form of this off-shell Coulomb singularity reveals that correc- 



tions of some per cent occur even a few W-decay widths above threshold. As explained 
above, the DPA becomes valid only several widths above threshold. Nevertheless, there 
exists an overlapping region in which the inclusion of the Coulomb singularity within the 
double-pole approach is reasonable. 

The Feynman graph relevant for the Coulomb singularity is diagram (d) of Fig. [l|. The 
non-factorizable corrections contain just the difference between the off-shell and the on- 
shell contributions of this diagram. Therefore, the difference between off-shell and on-shell 
Coulomb singularity is in principle included in A mm / , as defined in ( f4.4|) . The genuine form 
of A mm / in DPA, which is given by (|4.8|), does not contain the full effect of the Coulomb 



singularity, because in both Co functions of (f4.4j) the on-shell limit K± — > was taken 
under the assumption of a finite /?w- In order to include the correct difference between 
off-shell and on-shell Coulomb singularity in A mm /, the on-shell limit of the Cq functions 
of ( [4.4[ ) has to be taken for arbitrary f3yj. The full off-shell Coulomb singularity can be 
included by adding 



2M^-s 



2vri / /? + A M \ _ 2vri (K+ + K_ + s(3 w A M 



(3 \p + A M +pj Ay V 2fas 



(4.11) 



to the genuine DPA result ( |4.8|) for A mm /. The quantities /?, (3, and Am are defined in 
(A..1). After combination with the factorizable doubly-resonant corrections, all doubly- 
resonant corrections and the correct off-shell Coulomb singularity are included. The on- 
shell Coulomb singularity contained in the factorizable corrections is compensated by a cor- 
responding contribution in the non-factorizable ones. Note, however, that this subtracted 
on-shell Coulomb singularity appears as an artefact if the non-factorizable corrections are 
discussed separately. 
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4.4 Non-local cancellations 



In Ref. it was pointed out that the non-factorizable photonic corrections completely 
cancel in DPA if the phase-space integration over both invariant masses of the W bosons 
is performed. This cancellation is due to a symmetry of the non-factorizable corrections. 
The lowest-order cross section in DPA is symmetric with respect to the "reflections" 



(A* 



Ml) 

mi. 



+ iM w r w 
+ iM w r w 



{kl 
{kl 



Mir) 
Ml) 



iM w r w 
iM w r w 



-K* 
-K* . 



(4.12) 



Therefore, A can be symmetrized in 



-K* + or K_ 



-K*_ if the respective invariant 



mass is integrated out. For instance, if k_ is integrated out, A can be replaced by 



A + A 



~ 17T 



if_<->— if* 

s - 2Ml 



+ K + S 2 3Kw 



w / K_x w + K* 



K_xw — K + 
K_ K 



det(y) det(y) 



K--t) +ln ( l+ SM 



+ lnf 1 + ^(1 

V s 23 

i^_det(y 2 ) Kl det(K 2 ' 



det(y) 
^+det(F 3 



det(y) 
X + det(y 3 ') 



In hr w + 
In ( x w + 



det(y') 

S23 + -§24 



In (xy? + 



S23 

ln( I + ^,x w 



Sl3 + «23 



— In — x w 



■§23 



S23 + S 2 4 



+ 



det(y') 

2 g2 3M^(^ g2 4-^ 2 Sl3) 

det(y) 

2523^(^24-^X3 



Ml 



det(y) 
+ imaginary parts, 



lull - 
-lnfl + 



- 21nl 1 
-21nfl - 



K + S 23 + S24 

K- Ml 

K + 5 2 3 + S 2 A 

Kl Ml 



K-Xw 
Klx w 



(4.13) 



with z defined in ( |U.4j) . Note that this expression is considerably simpler than the full 
result for A. In particular, all dilogarithms have dropped out. 

Symmetrizing (|4.13|) also in K + <-> —Kl leads to further simplifications if the relations 
( |3.41| ) for the determinants are used. Under the assumptions that («X3 + S23) > -M^iw, 
(•§23 + ^24) > —Mlxw, and that k w is imaginary, the real part of the result vanishes. 
These assumptions are fulfilled on resonance, k\ = Ml; off resonance, there are regions 
in phase space where the assumptions are violated. The volume of those regions of phase 
space is suppressed by factors \k± — Ml\/Ml and thus negligible in DPA. Therefore we 
can use the above assumptions and find that A vanishes in DPA after averaging over the 
four points in the (k+, k 2 



plane that are related by the reflections (|4.12 ): 



A + A 


+ A 




f A 




K+^-Kl 


K-^-K*_ 





K±- 



-K1 



0. 



(4.14) 
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This explicitly confirms the results of Ref . M . In particular, the non-factorizable correc- 
tions vanish in the limit \k± — M^\ -C TwMw, i.e. for on-shell W bosons. 

The above considerations lead to the following simplified recipe for the non-factorizable 
corrections to single invariant-mass distributions, i.e. as long as at least one of the invariant 
masses k± is integrated out: the full factor A can be replaced according to 



A 



-|A + A 



-K* 



K 4 



-K* 



(4.15) 



where the first term on the r.h.s. is given in ( f4.13|) , and the second follows from the first 
by interchanging K + «-> K_ and S13 <-> S24. Note that no double counting is introduced, 
since the first (second) term does not contribute if k\ (k 2 _) is integrated out. In order 
to introduce the exact Coulomb singularity, one simply has to add the additional term of 
([in]) to fl£I5D and flCIgp . 



4.5 Non-factorizable corrections to related processes 

Since all non-factorizable corrections involving the initial e + e~ state cancel, the above 
results for the correction factor also apply to other W-pair production processes, such as 
77 — > WW — > 4 fermions and qq — > WW — > 4 fermions. 

The presented analytical results can also be carried over to Z-pair-mediated four- 
fermion production, e + e~ — > TIL — > 4 fermions. In this case, A ff / yields the complete non- 
factorizable correction, where all quantities such as M w and T w defined for the W bosons 
are to be replaced by the ones for the Z bosons. The fact that Q\ = Q2 and Q3 = Q4 
has two important consequences. Firstly, it implies the cancellation of mass singularities 
contained in A ff / when all contributions are summed as in ( |4.2|) . Secondly, it leads to the 
antisymmetry of <5 n f under each of the interchanges k\ <-> ki and k% ^ k^. It is interesting 

are directly applicable, since A m f and A n 



to note that ( f4.2|) with (|4.3|) are directly applicable, since A m f and A mm / cancel in the 
sum of (O). Therefore, a practical way to calculate the non-factorizable corrections to 
e + e~ — ► ZZ — > 4 fermions consists in taking A ff / + A mf / from (|4.9|) and ([4.10 ), and setting 



A, 



to zero. 



5 Numerical results 

If not stated otherwise, we used the parameters 

a' 1 = 137.0359895, M z = 91.187 GeV, M w = 80.22 GeV, T w = 2.08 GeV, (5.1) 

which coincide with those of Ref. M, for the numerical evaluation. 

In order to exclude errors, we have written two independent programs for the cor- 
rection factor Q4.2p and compared all building blocks numerically. These subroutines are 
implemented in the Monte Carlo generator EXCALIBUR [|1J as a correction factor to the 
three (doubly-resonant) W-pair-production signal diagrams. In all numerical results be- 
low, only these signal diagrams are included, and no phase-space cuts have been applied. 
We restrict our numerical analysis to the correction factor for purely leptonic final states, 
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e + e~ — > WW — > U£i + £'~ U£r . Results for other final states can easily be obtained from our 
formulae and computer programs. 

The results for the figures have been obtained from 50 million phase-space points using 
the histogram routine of EXCALIBUR with 40 bins for each figure. For each entry in the 
tables, 10 million phase-space points were generated. 



5.1 Comparison with existing results 

The non-factorizable photonic corrections have already been evaluated by two groups. As 
was noted in the introduction, the authors of Refs. [|, [L(J have not confirmed^ the results 
of Ref. ||. Therefore, we first compare our findings with the results of these two groups. 

Melnikov and Yakovlev || give the relative non-factorizable corrections only to the 
completely differential cross section for the process e + e~ — > WW — > iy e e + e~u e as a function 
of the invariant mass M + of the z/ e e + system for all other phase-space parameters fixed. 
All momenta of the final-state fermions lie in a plane, and the momenta k.2 and k3 of the 
final-state positron and electron point into opposite directions. The angle between the 
W~ boson and the positron is fixed to 6*w-e+ = 150°, and the CM energy is chosen as 
yfs = 180 GeV. The invariant mass of the e~u e system takes the values M_ = 78 and 
82 GeV. The other parameters are a = 1/137, M w = 80 GeV, and T w = 2.0 GeV. In 
Fig. ^ we show our results for the non-factorizable photonic corrections for this phase- 
space configuration. The intermediate-intermediate (mm') corrections agree reasonably 
with those of Ref. |§, but the other curves differ qualitatively and quantitatively^. We 
mention that Fig. || has been reproduced JHJ by the authors of Refs. || [L0| within the 



expected level of accuracy. While the individual contributions shown in Fig. |6] are at the 
level of 10% owing to mass singularities, the sum, which is free of mass singularities, is 
below 1.2%. 

Beenakker et al. || have evaluated the relative non-factorizable corrections to the 
distributions dcx/dM + dM_, da/dM + , and da/dM av , where M av = (M_ + M+)/2. Our 
corresponding results for the set of parameters given in (|5.1| ) are shown in Fig. [7| for the 
single invariant-mass distributions and in Table [l| for the double invariant-mass distribu- 
tion. The deviation between the distributions dcr/dM + and d<r/dM_ in Fig. ^|, which 
should be identical, gives an indication on the Monte Carlo error of our calculation. The 
single and double invariant- mass distributions agree very well with those of Ref. |J. The 
worst agreement is found for small invariant masses and amounts to 0.03%. In fact, the 
agreement is better than expected, in view of the fact that our results differ from those 
of Ref. P] by non-doubly resonant corrections. In the numerical evaluations of Ref. M 



the phase space and the Born matrix element are taken entirely on shell JT^]. Moreover 



3 While the result of Refs. j^, |lO) (as our result) for the complete non-factorizable correction is free of 
mass-singular logarithms, the result of Ref. || contains logarithms of ratios of fermion masses. However, 
the authors of Ref. || have informed us |l6| that the results of Ref. || and Refs. |l(J agree for equal 
fermion pairs in the final state. 

4 Although not stated in Ref. mass-singular parts have been dropped there in the numerical evalu- 
ation |l6| rendering a thorough comparison of the (mf') and (ff') parts impossible. Comparing the sum 
of all three contributions, i.e. the complete non-factorizable correction factor, our result differs from the 
sum read off from Ref. [B| . 
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Figure 6: Relative non-factorizable correction factor to the differential cross section for 
the phase-space configuration given in the text. The curves labelled mm', mf', and ff' cor- 
respond to the curves A, B, and C of Ref. |§, respectively. 
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Figure 7: Relative non-factorizable corrections to the single invariant-mass distributions 
da/dM± and da/dM av for the CM energy yfi = 184 GeV. 
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Table 1: Relative non-factorizable corrections in per cent to the double invariant-mass 
distribution dcx/dM + dM_ for the CM energy yfs = 184 GeV and various values of M± 
specified in terms of their distance from Mw in units of Tw, i.e. A± = (M± — Mw)/Fw- 



the scalar integrals are parametrized by scalar invariants different from ours, leading to 
differences of the order of \k± — M^\/M^. 

In Ref. 0, additionally, the decay-angular distribution dcr/dM_dM + dcos6'w+e+ has 
been considered, where #w+e+ is the decay angle between k + and k 2 in the laboratory 
system. Our results for this angular distribution are shown in Fig. ||. The cross section 
is small for cos#w+e+ ~ — 1, where the corrections are largest. Unfortunately the cor- 
responding figure in Ref. || is not correct flT7 |. The authors of Ref. have provided a 
corrected figure, which agrees reasonably well with Fig. but does not show the kinks 
in the curve for M± = 78 GeV. The kinks are due to a logarithmic Landau singularity 
in the 4-point functions. If one employs the on-shell parametrization of phase space, as 
in Ref. H], the Landau singularities appear at the boundary of phase space. Although no 
kinks appear in the physical phase space in this case, the Landau singularities still give 
rise to large corrections for cos#w+e+ ~ — 1- Since the kinks appear in a region where the 
cross section is small, they are not relevant for phenomenology. The issue of the kinks is 
further discussed in Sect. |5l 



5.2 Further results 

In Fig. H we show the non-factorizable corrections to the single invariant-mass distribution 
d<r/dM + for various CM energies. While the corrections reach up to 1.3% for y/s = 
172 GeV, they decrease with increasing energy and are less than 0.04% for ^fs = 300 GeV. 
Note that the shape of the corrections is exactly what is naively expected. If a photon 
is emitted in the final state, the invariant mass of the fermion pair is smaller than the 
invariant mass of the resonant W boson, which is given by the invariant mass of the 
fermion pair plus photon. Since we calculate the corrections as a function of the invariant 
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Figure 8: Relative non-factorizable corrections to the decay-angular distribution 
do-/dM_dM + dcos# w + e + for fixed values of the invariant masses M± and the CM energy 
v/i = 184 GeV. 
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Figure 9: Relative non-factorizable corrections to the single invariant-mass distribution 
dcr/dM + for various CM energies. 
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v/i/GeV 


172 


184 


192 


200 


300 


AM + /MeV 


-2.0 


-1.1 


-0.8 


-0.6 


-0.09 



Table 2: Shift of the maximum of the single invariant-mass distributions da/dM + induced 
by the non-factorizable corrections at various CM energies. 
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Figure 10: Relative non-factorizable corrections to the angular distribution 
d<r/dM_dM + dcos(/> for fixed values of the invariant masses M± and the CM energy 
y/s = 184 GeV. 



masses of the fermion pairs, the cross section tends to increase for small invariant masses 
and decrease for large invariant masses. 

The non-factorizable corrections distort the invariant-mass distribution and thus lead 
to a shift in the W-boson mass determined from the direct reconstruction of the de- 
cay products with respect to the actual W-boson mass. This shift can be estimated by 
the displacement of the maximum of the single-invariant-mass distribution caused by the 
corrections shown in Fig. |9|. To this end, we determine the slope of the corrections for 
M + = Mw, multiply this linearized correction to a simple Breit-Wigner factor, and de- 
termine the shift AM + of the maximum. The smallness of the correction allows us to 
evaluate AM + in linear approximation, leading to the simple formula 



AM, 



VdM, J 



r 



2 

W 



(5.2) 

M + =M W O 



Extracting the slope from our numerical results we obtain the mass shifts shown in Table |2|. 



In Figs. [10] and [TTJ we show the effect of the non-factorizable corrections on various 
angular distributions. Since the non-factorizable corrections are independent of the pro- 
duction angle of the W bosons, it suffices to consider distributions involving the angles 
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Figure 11: Relative non-factorizable corrections to the angular distribution 
da/dM_dM + dcos# e +e- for fixed values of the invariant masses M± and a CM energy 
y/s = 184 GeV. 



of the final-state fermions. We define all angles in the laboratory system, which is the 
CM system of the production process. The distribution over the angle <fi between the two 
planes spanned by the momenta of the two fermion pairs in which the W bosons decay, 
i.e. 

(ki x k 2 )(k 3 x k 4 ) 

cos0 = - — —, (5.3) 

|ki x k 2 ||k 3 x k 4 | 

is presented in Fig. [K]. The corrections are of the order of 1% or less. Like the (ft dis- 
tribution, the distribution over the angle between positron and electron ^e+e- (Fig. [TTJ) is 
symmetric under the interchange of M + and M_. As for the # e +w+ distribution (Fig. |]), 
the corrections reach several per cent in the region where the cross section is small. 



The distribution in the electron energy E e - is considered in Fig. 12. The corrections 



are typically of the order of 1%. Again the corrections become large where the cross section 
is small. 

In Sect. [4.3| we have introduced a correction term that includes the full off-shell Cou- 
lomb singularity. The results for the non-factorizable corrections with this improvement 



are compared with those of the pure DPA in Fig. |13[ For y 7 ? = 184 GeV the additional 



terms shift the non-factorizable corrections by up to 1.4% for dcr/dM av and by up to 
0.8% for dcr/dM + for small invariant masses, whereas for large invariant masses there is 
practically no effect. The difference originates essentially from the differences between 
1/(3 and l//3w i n (|4.11|) . For large invariant masses, the explicit logarithms in (^4.11|) are 
small, i.e. the Coulomb singularity correction is minuscule, and this difference practically 
makes no effect. For small invariant masses, the logarithms are approximately m and 
the difference causes the effect seen in Fig. [L3|. In Table ||] we show the non-factorizable 
corrections to the double invariant-mass distribution, as in Table [I], but now with the 
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Figure 12: Relative non-factorizable corrections to the electron-energy distribution 
da / dM_dM + dE e - for fixed values of the invariant masses M± and a CM energy yfs = 
184 GeV. 
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Figure 13: Relative non-factorizable corrections to the single invariant-mass distributions 
da/dM + and dcr/dM av for a CM energy y/s = 184 GeV with (lower curves) and without 
(upper curves) improved Coulomb-singularity treatment. 
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0.15 


0.10 


0.03 


-0.08 


-0.18 


-0.27 


-0.38 


1/2 


0.05 


-0.02 


-0.09 


-0.18 


-0.27 


-0.35 


-0.44 


1/4 


-0.04 


-0.13 


-0.19 


-0.27 


-0.35 


-0.41 


-0.49 


1 


-0.20 


-0.27 


-0.33 


-0.38 


-0.44 


-0.49 


-0.54 



Table 3: Same as in Table [T] but with improved Coulomb- singularity treatment. 



improved Coulomb-singularity treatment. We find a difference of up to half a per cent 
for small invariant masses but no effect for large ones. We mention that the difference 
between the entries in Tables [5] and [I] is directly given by the contribution ( 4.11| ) to A mm /, 
without any influence of the phase-space integration. 



5.3 Discussion of intrinsic ambiguities 

In the results presented so far, all scalar integrals were parametrized by s, s 2 3, S13, s 2 4, and 
k\ (parametrization 1). In DPA, however, the parameters of the scalar integrals are only 
fixed up to terms of order k\ — M^. We can for example parametrize the scalar integrals in 
terms of s, s 23 , Si 23 , s 234 , and k± (parametrization 2) instead. As a third parametrization, 
we fix all scalar invariants except for k± by their on-shell values, corresponding exactly to 
the approach of Ref. ||. The results of these three parametrizations differ by non-doubly- 
resonant corrections. 



The difference between parametrizations 1 and 2 is illustrated in Fig. |14] for the single 
invariant-mass distribution. The difference amounts to ~ 0.1%. Note that for an invariant 
mass M+ = 70 GeV we have a\M^ — k\\/M^ ~ 0.002 and would thus expect absolute 
changes in the non-factorizable corrections at this level. 

For the non-factorizable corrections to the angular distributions, uncertainties of the 
same order are to be expected. The only exceptions are the distributions over the decay 
angles #w+e+ and 6*w- e -. Let us explain this for 6*w+e+ in more detail: the non-factorizable 
correction contains the term 27riln[l + xw^w/( s i3 + ' s 23)], which can be evaluated by tak- 
ing (si3 + s 23 ) directly or (s 12 3 — M^) as input. This parametrization ambiguity can lead to 
larger uncertainties, because the above logarithm can become singular, and the location of 
this Landau singularity is shifted by the ambiguity. Since there is a one-to-one correspon- 
dence between S123 and #w+e+ for fixed s and k+, this logarithmic singularity is washed 
out if the angular integration over #w+e+ is performed, but appears as a kink structure in 
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Figure 14: Relative non-factorizable corrections to the single invariant-mass distribution 
da/dM + for the CM energies 172, 184, and 192 GeV and two different parametrizations. 
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Figure 15: Relative non-factorizable corrections to the decay-angular distribution 
dcT/dM_dM + dcos6' w + e + for fixed values of the invariant masses M± and a CM energy 
s/s = 184 GeV using three different parametrizations, as specified in the text. 
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the angular distribution over this angle. Figure 15 shows the non-factorizable corrections 
to this angular distribution for the three parametrizations. For M + = 78 GeV we still 
have a\k^_ — M^\/M^ ~ 0.0004. Apart from the regions where the Landau singularities 
appear, this is indeed of the order of the differences between the three parametrizations. 
When considering the Landau singularities, one should realize that the parametrization 
ambiguity of the locations of the singularities is not suppressed by a factor a, i.e. different 
parametrizations shift the locations at the level of \k+ — M^\/M^ ~ 0.05. However, the 
impact of the corresponding kinks on observables is again suppressed with M^|/M^ 
if the angles are integrated over, since the singularities can only appear near the bound- 
ary of phase space and disappear from phase space exactly on resonance. Since the cross 
section is small where the Landau singularity appears, the effect is phenomenologically 
irrelevant^. 



6 Conclusion 

A reasonable approach for the evaluation of the 0(a) corrections to W-pair-mediated 
four-fermion production consists in an expansion around the poles of the W propagators. 
For most of the corrections, the leading terms in this expansion, the doubly-resonant 
corrections, are sufficient. Combining the results for on-shell W-pair production and decay, 
one obtains the factorizable doubly-resonant corrections, which account for the largest 
part of the corrections. The remaining doubly-resonant corrections comprise the effects of 
virtual soft-photon exchange and real soft-photon interference where W-pair production 
and W decay are not independent. Those corrections are called non-factorizable. 

We have studied the non-factorizable doubly-resonant corrections to W-pair production 
analytically as well as numerically. By defining these corrections in such a way that the 
complete doubly-resonant corrections are obtained by adding the factorizable corrections, 
the non-factorizable corrections become manifestly gauge-independent. We have modified 
the double-pole approximation in such a way that the non-factorizable corrections also 
include the correct off-shell Coulomb singularity, which is the most important effect of the 
instability of the W bosons. 

It turns out that the actual form of the real non-factorizable corrections depends on the 
parametrization of phase space, i.e. on whether the photon momentum is included in the 
definition of the invariant masses of the resonant W bosons or not. This ambiguity is due 
to the fact that not only extremely soft photons but also photons with energies of the order 
of the W width contribute to the doubly-resonant corrections. For the parametrization 
where the invariant masses of the final-state fermion pairs are fixed when integrating over 
the photon momentum, we have listed all relevant analytical results. We have sketched 

Although it would be possible to remove all Landau singularities from the physical phase space in 
DPA by choosing an appropriate parametrization (e.g. parametrization 3), we have kept these structures, 
because we view this kind of ambiguity as unavoidable in DPA. Only a complete off-shell calculation in 
0(a) can resolve this uncertainty. We have checked that the Landau singularities are definitely present in 
the scalar integrals appearing in the exact off-shell calculation and we do not see a reason why these should 
drop out in the complete off-shell 0{a) calculation. In any case, the kink structures will be smeared out 
by the finite decay width of the W boson. 
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the calculation of the real bremsstrahlung integrals and the reduction of the real 5-point 
functions. The final analytical results are presented as transparent as possible so that 
they can be easily implemented in Monte Carlo generators. Our correction factor for the 
non-factorizable corrections agrees analytically with the one of Refs. || [10[] in double-pole 
approximation. 

We have numerically discussed various distributions and, in particular, compared to 
the results of Refs. H [| as far as possible. Our numerical results for the correction 
factor differ from those of Ref. ||, but agree with those of Ref. [jj at the level of the 
expected intrinsic uncertainty of the double-pole approach. These small deviations are 
due to different embeddings of the non-factorizable correction in the off-shell phase space, 
which are equivalent in double-pole approximation. 

For the adopted phase-space parametrization the non-factorizable doubly-resonant cor- 
rections are independent of the W-production angle, because all non-factorizable virtual 
and real soft-photonic effects that are related to the initial state cancel. For this reason, the 
presented correction factor for the non-factorizable doubly-resonant corrections is directly 
applicable to other reactions with the same final state, such as 77 — > WW — > 4fermions 
and qq — > WW — > 4fermions. The non-factorizable corrections to Z-pair-mediated four- 
fermion production, e + e _ — > ZZ — > 4fermions, can also be read off from our analytical 
results. 

The non-factorizable corrections are non-vanishing only for distributions in the in- 
variant masses of the W bosons and vanish after integration over both invariant masses. 
Therefore, they do not modify the usual genuine differential distributions that are used for 
the investigation of the triple gauge couplings. For invariant-mass distributions, the non- 
factorizable doubly-resonant corrections in the LEP2 energy region are small with respect 
to the experimental accuracy at LEP2, but should become relevant for possible future 
linear colliders with higher luminosity. At high energies the non-factorizable corrections 
are negligible. 
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Appendix 

A Useful definitions 

In the main text, we have already used the following short-hand expressions, 

Av = yi — — + ie, 
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x w 

p 
p 



Pw + l 



4M 2 



8 

y/Ks, kl k 2 



5i3-S24, S14S23) — ie, 



where 



A(x, y, z) = x 2 + y 2 + z 2 — 2xy — 2xz — 2yz 
The evaluation of one-loop n-point functions naturally leads to the usual dilogarithm 

dt 



(A.l) 
(A.2) 



r z rir 

U 2 (z) = - -ln(l-t), 
Jo t 

and its analytically continued form 



arc( 1 — z)\ < 7T, 



(A.3) 



£i 2 (x, y) = Li 2 (l - xy) + [ln(xy) - ln(x) - ln(y)] ln(l - xy), 

I arc(x)|, I arc(y)| < n. (A. 4) 

B Calculation of the real bremsstrahlung integrals 

In this appendix we describe a general method for evaluating the bremsstrahlung 3- and 



4-point integrals defined in ( 3.15 ). We make use of the generalized Feynman-parameter 



representation [IS 



1 f°° 

— = r(») y o dx,...^ 



B.l 



where the real variables oti > are arbitrary, but not all equal to zero simultaneously. The 
sum X)?=i N[xi must be non-zero over the entire integration domain. 

We first consider the 3-point integral Co- We only need the difference between the 
general IR-finite integral and the corresponding IR-divergent on-shell integral, 



C {p 1 ,p 2 ,0,m 1 ,m 2 ) - C (pi,p 2 ,X, ypj, Vpi 
■d 3 qj 1 



vrg I (2pig + pf - m\){2p 2 q + p 2 2 - m 2 2 ) (2p 1 q) (2p 2 q) 



;b.2) 



which is UV-finite and Lorentz- invariant. Extracting the signs cij of p^ via definition ( |3.29|) , 
and assuming for the moment Oi Re(p 2 — m 2 ) > 0, the Feynman-parameter representation 
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( |B.1|) can be applied to each term in the integrand of ( |B.2|) . The integration over d 3 q can 
be carried out and yields 



Cd(pi, P2, 0, mi, m 2 ) - C (pi,p 2 , A, Jp\ , Jp 2 2 



ln(A 2 ) + In 



;b.s) 



o"iO"2^ dxidx 2 5 (l -^aiX^j 



i=l 



2 In 



Y a i(P 2 i 



i=l 



2 x 2 

2(YPi X i 
=1 



Putting for instance a± = and a 2 = 1, and using 

1 1 



dx 



X + Z\ x + z 2 



ln(l + zx) = Ci 2 {zi, z) - Ci 2 (z 2 , z) 



(B.4) 



where | arcf^i^, z, z zi j2 )\ < it, the remaining one- dimensional integration and the analyti- 
cal continuations to arbitrary complex (p 2 — mf) are straightforward. 
Next we consider the IR-fmite 4-point function 



V (p 1 ,p 2 ,p 3 ,0,m 1 ,m 2 ,m 3 ) 
"d 3 q 



irq (2p 1 q +p{- m\)(2p 2 q + p 2 2 - m 2 2 )(2p 3 q + pj - m| 



For <Tj Re(pf — m 2 ) > we can proceed as above and find 



(B.5) 



<3o=M 



A)(Pi,P2,P3,0,rai,m 2 ,m 3 ) = a x a 2 a 3 / dxi<lx 2 dx. 



i=l 



r 3 



i=l 



53 Pt 

i=l 



(B.6) 

Again, the two-dimensional integration over the Feynman parameters and the analytical 
continuations in (p 2 — m 2 ) are straightforward. 

Finally, we inspect the IR-divergent 4-point function 



£>o(pi,P2,P3, A, \fpj,m 2 , y^) = J 



1 



d 3 q 

vrg {2piq)(2p 2 q + p 2 2 - ml){2p 3 q) 



(Jo- 



V c i 2 + A2 



;b.7) 



Instead of applying the Feynman-parameter representation directly, it is more convenient 
to extract the IR singularity by subtracting and adding an IR-divergent 3-point function 
containing the same IR structure. Since the difference between the 4- and 3-point functions 
is IR-finite, we can regularize the IR divergence in these two integrals in a more convenient 
way. Therefore, we write 



^o(Pi,P2,P3, A, Jpf,m 2 , slpi) = lim 



m 2 — >p 2 



A) (Pi,P2,P3,0,m,m 2 , Jp% 



—5— — 2 C o(pi,P3,0,m, JpJ 
pi -mi v v 



1 



2 Co[Pi,P3,X, y Pi, VPs 



(B.8) 
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i.e. we regularize the IR divergence in the 4-point function by the off-shellness pf — m 2 ^ 0. 
Both the 4-point function as well as the difference of 3-point functions can be treated 
as above, yielding straightforward Feynman-parameter integrals according to ( |B.3| ) and 
( B.6|) , respectively. The limit m 2 —>■ p\ can easily be taken after the integrations have 
been performed. 



C Explicit results for scalar integrals 

C.l Loop integrals in double-pole approximation 

In the (mm') correction the following combination of 3-point functions appears in the strict 
DPA: 



C (k + , -k-, 0, M, M) - [C (k+, —k-,X, M w , Af w ) 

^-S^- Ci 2 (^,x w ^j + Ci 2 (j^,x^fj +Li 2 (l-a;w) + 7r 2 + ln 2 (-xw) 
ln(xw) — 27riln(l — #w) j- 



+ 2 In 



K 4 



(C.l) 



In order to include the full off-shell Coulomb singularity, one has to add 



27Ti / {3 + Am-P 
sfi n \f3 + A M +P 



\ _ 2vri_ ( K + + K^ + s(3 w A M \ 



to the r.h.s. of Q. 

For the (mf') and (ff') corrections we need the following 4-point functions: 



(C.2) 
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(C.3) 
(C.4) 

(C.5) 



-D (2) = Do(-h, k+, k 2 , A, m 3 , M, m 2 ) 
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(C.6) 
(C.7) 
(C.8) 



C.2 Bremsstrahlung integrals in double-pole approximation 

In the (mm') interference corrections the following combination of 3-point functions ap- 
pears: 



kl=Ml 
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(C.9) 



For the (mf') and (ff') interference corrections the following 4-point functions are re- 
quired: 
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K-^-K* R*_S 23 



iK*_m 2 



A(s 2 3 + S24) 



(C.11) 



(C.12) 



(C.13) 



P (4) = ^ofe fc-, 0, m 3 , M*, M) 
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-A)(4) 



K- 



+ 



2vri 



lnfl + XwM w 



■ 



Sl3 + S23J - 

(C.14) 



We note that the logarithmic terms on the r.h.s. of (|C.10| ) yield purely imaginary contri- 
butions if (S13 + S23) > -M^,x w , (s 2 3 + s 24 ) > —MyjXwi an d is imaginary. These 
conditions are fulfilled on resonance. Off resonance, this is no longer true near the bound- 
ary of phase space. But since these conditions are only violated in a fraction of phase 
space of order \k± — M^\/M^, which is irrelevant in DPA (compare the discussion of the 
relevance of the Landau singularities in Sect. |5.3| ), it would even be allowed to replace 
A)(0), 



which is real, by — Re{D (0)} , as it was done in Refs. [BL ITU 
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